Pairing of bosons in the condensed state of the boson-fermion model 
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A two component model of negative (7 centers coupled with the Fermi sea of itinerant fermions is 
discussed in connection with high-temperature superconductivity of cuprates, and superfluidity of 
atomic fermions. We examine the phase transition and the condensed state of this boson-fermion 
model (BFM) beyond the ordinary mean-field approximation in two and three dimensions. No 
pairing of fermions and no condensation are found in two-dimensions for any symmetry of the order 
parameter. The expansion in the strength of the order parameter near the transition yields no 
linear homogeneous term in the Ginzburg-Landau-Gor'kov equation and a zero upper critical field 
in anj/-dimensional BFM, which indicates that previous mean-field discussions of the model are 
flawed. Normal and anomalous Green's functions are obtained diagrammatically and analytically 
in the condensed state of a simplest version of 3D BFM. A pairing of bosons analogous to the 
Cooper pairing of fermions is found. There are three coupled condensates in the model, described 
by the off-diagonal single-particle boson, pair-fermion and pair-boson fields. These results negate 
the common wisdom that the boson-fermion model is adequately described by the BCS theory at 
weak coupling. 

PACS numbers: PACS: 71.20.-z,74.20.Mn, 74.20.Rp, 74.25.Dw 



I. INTRODUCTION 



The experimental 0, 0, IE II IE H and theoretical 
0, IE EI El evidence for an exceptionally strong electron- 
phonon interaction in novel superconductors is now so 
overwhelming that even some advocates of the non- 
phononic mechanisms hdj admit the fact. A few authors 
(see, for example 0, 0, Q, El) explored a view that 
the extension of the BCS theory towards the strong in- 
teraction between electrons and ion vibrations describes 
the phenomenon. In this regime, pairing takes place in 
real space due to a polaron collapse of the Fermi en- 
ergy or due to a low density of carriers. At first 
sight, bipolarons have a mass too large to be mobile. 
Indeed, Anderson introduced small bipolarons as en- 
tirely localised objects explaining some unusual proper- 
ties of chalcogenide glasses. However, it has been shown 
more recently that the inclusion of the on-site Coulomb 
repulsion leads to the favoured binding of intersite carri- 
ers ^E EiJl ■ The intersite bipolarons can tunnel with an 
effective mass of about 10 electron masses fl5|. fl^ . flfll l20j . 
and account for a high critical temperature |21| . 

Soon after Anderson an d Street and Mott 
introduced localized pairs in amorphous semiconduc- 
tors, a two component model of negative U centers cou- 
pled with the Fermi sea of itinerant fermions was em- 
ployed to study superconductivity in disordered metal- 
semiconductor alloys |U 0] ■ When the attractive po- 
tential U is large, the model is reduced to localized hard- 
core bosons spontaneously decaying into itinerant elec- 
trons and vice versa, different from a non-converting mix- 
ture of mobile charged bosons and fermions [25, 26] . This 
boson-fermion model (BFM) was applied more generally 
to describe pairing electron processes with localization- 
delocalization |27 |. an d a linear resistivity in the normal 
state of cuprates j23|. The model attracted more atten- 
tion in connection with high-temperature superconduc- 



tors [1C1|29|,|^ |^|321|MIMI^IMIM|38|, as an alter- 
native to the bipolaron model of cuprates with intrinci- 
cally mobile bipolarons 39]. In particular, Refs. [3ril3^ ] 
claimed that 2D BFM with immobile hard-core bosons 
is capable to reproduce some physical properties and the 
phase diagram of cuprates. BFM has been also adopted 
for a description of superfluidity of atomic fermions scat- 
tered into bound (molecular) states |40|. 

Most studies of BFM below its transition into a low- 
temperature condensed phase applied a mean-field ap- 
proximation (MFA) , replacing zero- momentum boson op- 
erators by c-numbers and neglecting the boson self- 
energy in the density sum rule [^E I^E HE HE HE HE 
l38l yjj ■ When the bare boson energy is well above the 
chemical potential, the BCS ground state was found with 
bosons being only virtually excited |2E HE] • MFA led to 
a conclusion that BFM exhibits features compatible with 
BCS characteristics 32], and describes a crossover from 
the BCS-like to local pair behaviour [3^ . The transition 
was found more mean-field-like than the usual Bose con- 
densation, i.e. characterized by a relatively small value 
of the fluctuation parameter Gi [33|. 

At the same time our previous study of BFM HE] be- 
yond MFA revealed a crucial effect of the boson self- 
energy on the normal state boson s pec tral function and 
the transition temperature T c . Ref.|34| proved that the 
Cooper pairing of fermions via virtual bosonic states 
is impossible in any-dimensional BFM. It occurs only 
simultaneously with the Bose-Einstein condensation of 
real bosons. The origin of this simultaneous condensa- 
tion lies in a softening of the boson mode at T = T c 
caused by its hybridization with fermions. The energy 
of zero-momentum bosons is renormalized down to zero 
at T = T c , no matter how weak the boson-fermion cou- 
pling and how large the bare boson energy are [34j. One 
can also expect that the boson self-energy should quali- 
tatively modify the phase transition and the whole con- 
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densed phase of BFM below T c . 

In this paper the phase transition and the condensed 
state of BFM are examined beyond the ordinary mean- 
field approximation in two and three dimensions. It is 
shown that T c = K in the two-dimensional model, even 
in the absence of any Coulomb repulsion, and the phase 
transition is never a BCS-like second-order phase tran- 
sition even in 3D BFM because of the complete boson 
softening. A closed set of equations for fermion and bo- 
son Green's functions (GFs) is derived taking into ac- 
count the self-energy effects in the condensed state of 3D 
BFM. There exist a boson pair condensate along with 
the fermion Cooper pair and the single-particle boson 
condensate in the model. Remarkably, the Gor'kov ex- 
pansion j^H of GFs in the strength of the order parameter 
yields a zero linear term at any temperature below T c , 
and zero upper critical field. 



II. NO COOPER PAIRING AND 
CONDENSATION IN 2D BFM 

2D BFM is defined by the Hamiltonian, 

h = E ac k , CT c k , ff + £ $> q fo q + (i) 

k,o-=T4 q 

gN~ 1/2 Y^ (0kfr q C_ k +qAT C k+q/2,l +H.C.) , 

q,k 

where £k = ~2t(cosk x + cosk v ) — fJ, is the 2D energy 
spectrum of fermions, E = — 2fi is the bare boson 
energy with respect to their chemical potential 2/i, g is 
the magnitude of the anisotropic hybridization interac- 
tion, 0k = 0_k is the anisotropy factor, and N is the 
number of cells. Here and further I take fi=c=/cB = l 
and the lattice constant a = 1. Ref. [3(| argued that 
'superconductivity is induced in this model from the 
anisotropic charge-exchange interaction (g</>k) between 
the conduction-band fermions and the immobile hard- 
core bosons', and 'the on-site Coulomb repulsion com- 
petes with this pairing' reducing the critical temperature 
T c less than by 25%. Also it has been argued [33, that 
the calculated upper critical field of the model fits well 
the experimental results in cuprates. 

Here I show that T c = K in the two-dimensional 
model, Eq.(l), even in the absence of any Coulomb re- 
pulsion, and the mean-field approximation is meaningless 
for any-dimensional BFM because of the complete boson 
softening. 

Replacing boson operators by c-numbers for q = in 
Eq.(l) one obtains a linearised BCS-like equation for the 
fermion order-parameter (the gap function) A k , 



»k \ ^ 



k' 



A k , tanh(&'/2T) 



(2) 



with the coupling constant g 2 = g 2 (l — 2n B ), renor- 
malized by the hard-core effects. Using a two-particle 



fermion vertex part in the Cooper channel one can prove 
that this equation is perfectly correct even beyond the 
conventional non-crossing approximation [34[ . The prob- 
lem with MEA does not stem from this BSC-like equa- 
tion, but from an incorrect definition of the bare boson 
energy with respect to the chemical potential, Eq(T). 
This energy is determined by the atomic density of 
bosons (n B ) as (Eq.(9) in Ref. 



tanh 



Eq 

2T 



1 -2n 



B 



(3) 



While Eq.(2) is correct, Eq.(3) is incorrect because the 
boson self-energy £&(q, Q) due to the same hybridization 
interaction is missing. At first sight |3g the self-energy 
is small in comparison to the kinetic energy of fermions, 
if g is small. However £(,(0,0) diverges logarithmically 
at zero temperature |34j . no matter how week the in- 
teraction is. Therefore it should be kept in the density 
sum- rule, Eq.(3). Introducing the boson Green's function 



1 - 2n J 



ifl - E a - E 6 (q, fi) 
one must replace incorrect Eq.(3) by 



(4) 



T 

N 



(5) 



q,n 



where r = +0, and Q = 2nTn (n = 0, ±1, ±2...). 

The divergent (cooperon) contribution to S&(q, f2) is 
given by Fig.l 0, 



S 6 (q,0) 



,v2 



2N 



(0) 



tanh[£ k - q/2 /(2T)] + tanh[£ k+q /2/(ZT)] 



Ck-q/2 + Ck+q/2 ~ ^ 



so that one obtains 



E b (q,0) = S b (0,0) + 



2M* 



0(q 4 



(7) 



for small q and any anisotropy factor compatible with 
the point-group symmetry of the cuprates. Here M* is 
the boson mass, calculated analytically in Ref.[34| for 
the isotropic exchange interaction and parabolic fermion 
band dispersion (see also Ref.|3^]). The BCS-like equa- 
tion (2) has a nontrivial solution for Ak at T = T c , if 



E a = -E 6 (0,0). 



(8) 



Substituting Eq.(7) and Eq.(8) into the sum-rule, Eq.(5), 
one obtains a logarithmically divergent integral with re- 
spect to q, and 



con fL = Q _ 



Jo d lll 



(9) 
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of the self-energy) as |34 




FIG. 1: Diagram for the boson self-energy. Solid lines are the 
fermion GFs. Vertex (dot) corresponds to the hybridization 
interaction. 



The devastating result, Eq.(9) is a direct consequence of 
the well-known theorem, which states that BEC is im- 
possible in 2D. 

One may erroneously believe that MFA results [3(| |3tJ 
can be still applied in three-dimensions, where BEC is 
possible. However, increasing dimensionality does not 
make MFA a meaningful approximation for the boson- 
fermion model. This approximation leads to a naive 
conclusion that a BCS-like superconducting state oc- 
curs below the critical temperature T c ~ /iexp (~Eq/z c ) 
via fermion pairs being virtually excited into unoccupied 
bosonic states [2<||3(3- Here z c = g 2 N(0) and 2V(0) is the 
density of states (DOS) in the fermionic band near the 
Fermi level fi. However, the Cooper pairing of fermions is 
impossible via virtual unoccupied bosonic states also in 
3D BFM. Indeed, Eqs.(2,8) do not depend on the dimen- 
sionality, so that the analytical continuation of Eq.(4) 
to real frequencies ui yields the partial boson DOS as 
p(u) = (1 — 2riB)S(uj) at T = T c and q = in any- 
dimensional BFM. The Cooper pairing may occur only 
simultaneously with the Bose-Einstein condensation of 
real bosons in the exact theory of 3D BFM [34|. The 
origin of the simultaneous condensation of the fermionic 
and bosonic fields in 3D BFM lies in the softening of the 
boson mode at T = T c caused by its hybridization with 
fermions. 

Taking into account the boson damping and dispersion 
shows that the boson spectrum significantly changes for 
all momenta. Continuing the self-energy, Eq.(6) to real 
frequencies yields the damping (i.e. the imaginary part 



7(q,w) 



4q£ 



In 



cosh(g£ + u;/(4T c )) 



cosh(-<z£ + w/(4T c )) 



(10) 



where £ = vf/(4T c ) is a coherence length, and vf is 
the Fermi velocity. The damping is significant when 
q£ << I. In this region j(q,u>) = W7rz c /(8T C ) is com- 
parable or even larger than the boson energy lu. Hence 
bosons look like overdamped diffusive modes, rather than 
quasiparticles in the long- wave limit [3 HH > contrary to 
the erroneous conclusion of Ref . [3ll| . that there is 'the 
onset of coherent free-particle-like motion of the bosons' 
in this limit. Only outside the long- wave region, the 
damping becomes small. Indeed, using Eq.(10) one ob- 
tains 7(q, cj) = uiTz c /(2qvF) « w, so that bosons at 
q >> l/£ are well defined quasiparticles with a logarith- 
mic dispersion, u)(q) = z c \n(q£) 34]. Hence the boson 
energy disperses over the whole energy interval from zero 
up to Eq. 

The main mathematical problem with MFA in 3D also 
stems from the density sum rule, Eq.(5) which determines 
the chemical potential of the system and consequently 
the bare boson energy Eq(T) as a function of tempera- 
ture. In the framework of MFA one takes the bare boson 
energy in Eq.(2) as a temperature independent param- 
eter, E Q = g 2 N(0)\n(n/T c ) 33], or determines it from 
the conservation of the total number of particles neglect- 
ing the boson self-energy, Eq.(3) [M IS H HJ . Then 
Eq.(2) looks like the conventional linearized Ginzburg- 
Landau-Gor'kov equation |4l| with a negative coefficient 
a oc T — T c at T < T c in the linear term. Then 
one concludes that the phase transition is almost the 
conventional BCS-like transition, at least at Eq 3> T c 
[291130 . 33]. These findings are mathematically and phys- 
ically erroneous. Indeed, the term of the sum in Eq.(5) 
with Q n — is given by the integral 



T 



f/q 



1 



27r 3 £ + S b (q,0)' 



(11) 



The integral converges, if and only if Eq — Ef,(0, 0). In 
fact, 



£ + £ b (0,0) = 



(12) 



is strictly zero in the Bose-condensed state, because 
Hb = —[Eo + £b(0, 0)] corresponds to the boson chemical 
potential relative to the lower edge of the boson energy 
spectrum. More generally, /!/, = corresponds to the 
appearance of the Bogoliubov-Goldstone mode due to a 
broken symmetry below T c . This exact result makes the 
BSC equation (2) simply an identity [34| with a{T) = 
at any temperature below T c . On the other hand, MFA 
violates the density sum-rule, predicting the wrong neg- 
ative ot(T) below T c . 

Since a(T) — 0, one may expect that the conventional 
upper critical field, H C 2(T) is zero in BFM. To determine 
H c2 (T) and explore the condensed phase of 3D BFM, one 
can apply the Gor'kov formalism 41], as described below. 



4 



III. NORMAL AND ANOMALOUS GREEN'S 
FUNCTIONS OF 3D BFM: PAIRING OF BOSONS 

Let us now explore a simplified version of 3D BFM 
in an external magnetic field B = V x A neglecting the 
hard-core effects, 



H.c. 



(13) 



where ijj s (r) and 4>(r) are fermionic and bosonic fields, 
s =T, | is the spin, h(r) = — [V + ieA(r)] /(2m) — fx is 
the fermion kinetic energy operator. Here the volume of 
the system is taken as V — 1. 

The Matsubara field operators, Q = 
exp(i?r)Q(r)exp(~i?r),g = exp(iJr)gt(r) exp(-iJr) 
(Q = ip s ,<t>) evolve with the imaginary time 
-1/T r 1/T as 



Ot 

dr 
d<f>(r,r) 
8t 



h(r)ih (r, r) + <?0(r, t)^ 1 (r,r), (14) 
r(r)^(r,T)-#(r,r)^ T (r,r),(15) 
^0(r,r)+^i(r,r)V^(r,T). (16) 



The theory of the condensed state can be formu- 
lated with the normal and anomalous fermion GFs 
ELafciV) = -(T^ s (r,r)^(r',0)>, ^+(r,r',r) = 
{T r tp i(r, r)tp-i (r' , 0)}, respectively, where the operation T T 
performs the time ordering. Fermionic and bosonic fields 
condense simalteneously [34j. Following Bogoliubov ^2] 
the bosonic condensate is described by separating a large 
matrix element <^o( r ) m 4>{ v t t ) as a number, while the 
remaining part </>(r, r) describes a supracondensate field, 
(j)(r, t) = 0o ( r ) + <p(v, t). Then using Eq.( 16) one obtains 



#o(r) = A(r) 



E, 



-jF(r,r,0- 



(17) 



where !F{t,t',t) = (T t i/j i(r, r)-0f (r' , 0)}. The equations 
for GFs are obtained by using Eqs. (14-16) and the di- 
agrammatic technique l43l in the framework of the non- 
crossing approximation |4J|j as shown in Fig. 2 and Fig. 3. 

An important novel feature of BFM is a pairing of 
supracondensate bosons, caused by their hybridization 
with the fermionic condensate, as follows from the last 
diagram in Fig. 3. Hence, one has to introduce an 
anomalous supracondensate boson GF, S + (r, r',r) = 
{T T (f>(r,T)<f>(r' , 0)) along with the normal boson GF, 

P(r,r',r) = -(T T <f>(r, r)0(r',O)). The diagrams, Fig.2 
and Fig. 3, are transformed into analytical equations for 
the time Fourier-components of the fermion GFs with the 
Matsubara frequencies oj — 7rT(2n+l) (n = 0, ±1, ±2, ...) 

as 

[iw - h(T)]g u (r, r') = S(r - r') - A(r)jT+(r, r') 



» > = > + => ^ 

G 



+ 



> — — — ^- + > — — y% r - — 



-<■ = 

F + 



< > 



< frS J^C ^ + < ^ ^ ^ 



FIG. 2: Diagrams for the normal and anomalous fermion GFs. 
Zig-zag arrows represent the single-particle Bose condensate 
<f>o, dotted lines are the boson GFs. 



~ 9 2T Y.J ^_^(x,r)X» w _ w ,(r,x)g a; (x,rO 

UJ 1 

- 9 2 TY, J rfx^+(r,x)^+^(r,x)^+(x,r'), (18) 

U)' 

Ho;-^(r)]^(r ) r') = A*(r)^(r,rO 

- 9 2 TJ2 I dx£ w ,(r,x)P^(x,r)jF+(x,r') 

+ 9 2 TJ2J dx^,(r,x)^_ u ,(r,x)^(x,r') 



and, 



(iSl - E )V n {r,r') = 8{v - r') 
~ 9 2 TJ2J dxa w ,(r,x)^_^(r,x)D n (x,r') 

- 9 2 TY, J dx^(r,x)^j W (r,x)B+(x,r'),(19) 

U)' 

(-itl- E Q )B+ (r,r') = 

5 2 T^ J dx^,(r,x)^V^( r ' x )^(x,r') 
- ^ T E / dx^(x,r)^_o(x,r)B+(x,r'). 
for the boson GFs with B(t,t',t) = (T T ^(r,r)0(r', 0)). 

IV. GOR'KOV EXPANSION 

These equations can be formally solved in the homo- 
geneous case without the external field, A = 0. Trans- 
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s- ■> = > ■=>- + >■■■■ 

D 



o 



■< - = « <• 



o 



+ < <# - 



FIG. 3: Diagrams for the supracondensate boson GFs. The 
Cooper-pairing of fermions leads to the Cooper-pair-like bo- 
son condensate, described by the boson anomalous GF, B + . 

FIG. 4: The fermion self-energy 



forming into the momentum space yields GFs' time-space 
Fourier components as 



and 



0(k,w) 
J*-(k,w) 

D(q,w) 
B+(q,w) 



|i<D-al 2 + |A(k,o;)P 
A*(k,w) 



|^-£ | 2 + |r(q,ft)| 2 ' 
r»( q ,n) 
|zft-£ | 2 + |r(q,f>)| 2 ' 



(20) 
(21) 

(22) 
(23) 



where uj = uj + iS/(k, w), f2 = Q + iSfe(q, O), and 
£k = k 2 /(2m) —fx. The fermionic order parameter, renor- 
malised with respect to the mean-field A due to the for- 
mation of the boson-pair condensate, is given by 

A(k,c)= A + . 9 2 T]T J ^ + (k-q,0 J ')B(q,W + w'), 

(24) 

and the boson-pair order parameter, generated by the 
hybridization with the fermion Cooper pairs, is 



have about the same magnitude, as the fermion, Fig. 4, 
and boson, Fig.l, self-energies, 



E/(k,a;) 



-9 



rfq 

2^3 



^(q-k, -a/)£>(q,u;-a/), 
(26) 

E b (q,0) = -.g 2 T^ / 0£(k, o/)S(q - k, ft - a/), 

(27) 

respectively. 

On the other hand, when the temperature is close to 
T c (i.e. T c — T <C T c ), the boson pair condensate is 
weak compared with the single-particle boson and the 
Cooper pair condensates. In this temperature range T, 
Eq.(25) is of the second order in A, r oc A 2 , so that the 
anomalous boson GF can be neglected, since A is small. 
The fermion self-energy, Eq.(26) is a regular function of 
uj and k, so that it can be absorbed in the renormalizcd 
fermion band dispersion. Then the fermion normal and 
anomalous GFs, Eqs. (20,21) look like the familiar GFs of 
the BCS theory, and one can apply the Gor'kov expansion 
|41| in powers of A(r) to describe the condensed phase 
of BFM in the magnetic field near the transition. Using 
Eq.(18) one obtains to the terms linear in A 



r(q,fi) =9 2 TJ2 I ^(k,a/)^(q-k,Q-a/). ( 25 ) A*(r) = ^-Tj^ J dx(?i"J„ (x, r)A*(x)^ (x, r). (28) 



Hence, there are three coupled condensates in the model 
described by the off-diagonal fields g<j>o, A, and T, rather 
than two, as in MFA. At low temperatures all of them 



The spatial variations of the vector potential are small 
near the transition. If A(r) varies slowly, the normal 

state GF, Qw t \v, r') differs from the zero-field normal 
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state GF, ^ 0) (r - r') only by a phase [4l| gi n) (r,r') = 
exp[— ieA(r) • (r — i"')]^ -* (r — r'). Expanding all quan- 
tities near the point x = r in Eq.(28) up to the second 
order in x — r inclusive, one obtains the linearised equa- 
tion for the fermionic order parameter as 

7[V - 2ieA(r)] 2 A(r) = aA(r), (29) 

where 

a = l + ^M„i_ffl9) In M (30) 
and 7 « 7C(3)w|, 5 2 iV(0)/(487r 2 T 2 £;o)- 

V. CONCLUSION 

The coefficient a(T) disappears in Eq.(29), since Eq = 
— £(,(0,0) at and below T c , Eq.(12). It means that the 
phase transition is never a BCS-like second-order phase 
transition even at large Eq and small g. In fact, the tran- 
sition is driven by the Bose-Einstein condensation of real 
bosons with q = 0, which occur due to the complete soft- 
ening of their spectrum at T c . Remarkably, the conven- 
tional upper critical field, determined as the field, where 
a non-trivial solution of the linearised Gor'kov equation 
(29) occurs, is zero in BFM, H c2 (T) = 0. It is not a finite 
H c2 (T) found in Ref. JH using MFA. 

This qualitative failure of MFA might be rather un- 
expected, if one believes that bosons in Eq.(l) play the 



same role as phonons in the BCS superconductor. This 
is not the case for two reasons. The first one is the den- 
sity sum- rule, Eq.(5), for bosons which is not applied to 
phonons. The second being that the boson self-energy 
is given by the divergent (at T = 0) Cooperon diagram, 
while the self-energy of phonons is finite at small cou- 
pling. 

In the homogeneous case A(T) should be determined 
from Eq.(5) rather than from the BCS-like equation (2), 
which is actually the identity. To get an insight into 
the magnetic properties of the condensed phase one has 
to solve Eqs.(18,19) and Eq.(5) keeping the non-linear 
terms. Even at temperatures well below T c the con- 
densed state is fundamentally different from the MFA 
ground state, because of the pairing of bosons. The lat- 
ter is similar to the Cooper-like pairing of supraconden- 
sate 4 H e atoms (4^ , proposed as an explanation of the 
small density of the single-particle Bose condensate in 
superfluid Helium-4. The pair-boson condensate should 
significantly modify the thermodynamic properties of the 
condensed BFM compared with the MFA predictions. 
The common wisdom that at weak coupling the boson- 
fermion model is adequately described by the BCS the- 
ory, is therefore negated by our theory. 

I highly appreciate enlightening discussions with A.F. 
Andreev, L.P. Gor'kov, V.V. Kabanov, A.I. Larkin, A. P. 
Levanyuk, R. Micnas, and S. Robaszkiewicz, and support 
by the Leverhulme Trust (UK) via Grant F/00261/H. 
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